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1. Introduction 

Eugene Wigner introduced random matrix models about fifty years ago into nuclear 
pliysics [I|. Tfie scope of applications has increased over the years [2l [3] including 
fields such as molecular and atomic physics, mesoscopics and field theory. More 
recently random matrix theory has started to be used in quantum information theory 
[H El El [71 [8|. For an introduction to such applications see |9]. There the concept of 
individual qubits and their interactions becomes important. This implies that we enter 
the field of two-body random ensembles (TBRE) [TOl[lT], i.e. ensembles of Hamiltonians 
of n-body systems interacting by two-body forces. While such ensembles have received 
considerable attention, it was first focussed on fermions and later also included bosons. 
Yet in quantum information theory the qubits are taken to be distinguishable, and 
indeed the same holds for spintronics. Interest in both fields has sharply increased 
recently [121 [I3]. 

It is thus very pertinent to formulate and investigate TBRE's for distinguishable 
qubits. As random matrix ensembles are mainly determined by their symmetry 
properties this ensemble will be very different from other TBRE's. In particular, as 
the particles are distinguishable, their interaction can vary from particle pair to particle 
pair and can indeed be randomly distributed, thus introducing an entirely new aspect. 
This has the consequence that the topology according to which spins or qubits are 
distributed or interact will be important. Thus chains, trees and crystals of particles 
with nearest, second nearest and up to kth order interaction can be represented. 

As mentioned above, random matrix ensembles are usually basically defined by the 
invariance group of their measure and, if that is not enough, some minimal information 
conditions [HI Ej or independence condition [15] . Note that we deal with a symmetry of 
the ensemble, rather than with a symmetry of individual systems. The two concepts are 
to some degree complementary, and the former has also been called structural invariance. 
We propose an adequate definition for such ensembles in a very general framework 
in terms of independent Gaussian distributed variables. We then give an alternate 
representation in terms of the invariance group and variables that determine the orbits 
of the Hamiltonian on the ensemble under the action of the group. 

In order to show the relevance of the new ensemble we address the simplest possible 
topology, namely the chain with nearest neighbour interactions. For this system we focus 
on the ensemble averaged structure of the ground state and demonstrate the existence of 
an unusual quantum phase transition [16], which is triggered by breaking of time-reversal 
invariance (TRI). 

Entanglement, a key resource of quantum many-body systems in terms of quantum 
information, is to large extent related to quantum correlations, localization properties 
and quantum chaos. Entanglement has also been used as a property, alternative to 
long-range order in spatial correlation functions, to describe systems undergoing a 
quantum phase transition [17J. In one- dimensional systems such as quantum spin 
chains, it was shown [18] that the entanglement entropy of the ground state typically 
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saturates or diverges logarithmically with size when approaching the thermodynamic 
limit. Furthermore, it has been shown that logarithmic divergence implies quantum 
criticality. 

Interesting results emerge when a spatially homogeneous spin model is replaced by 
its disordered counterpart, where the spin interactions are taken at random. In this 
case there is often no physical justification why random interactions should still obey 
specific restricted forms such as Ising or Heisenberg interactions. In this context, we 
argue, it is more natural to use two-spin random ensembles (TSRE) for distinguishable 
particles, specifically choosing quantum spins 1/2, though these ensembles can readily 
be generalized to arbitrary spin. By construction these ensembles, as given in section 
[21 are invariant with respect to arbitrary local rotations, which we may view as gauge 
transformations. Another physical motivation for the definition of such ensembles is 
the coupling among arbitrary and perhaps mutually independent two level quantum 
systems which may come from completely different physical contexts such as e.g. two- 
level atoms, Josephson junctions and photons. 

In section |3] we concentrate on one-dimensional systems or spin chains and present 
results of numerical calculations, mainly based on density matrix renormalization group 
(DMRG) [19], in which we investigate entanglement and correlation properties of the 
ground state, averaged over an ensemble, and the average spectral gap to the first excited 
state as well as its fiuctuations. If we include the interaction with an external random 
magnetic field, and hence TRI is broken, we find fast decay of correlations, saturation 
of entanglement entropy, and power law decay g ~ A^"*^-^ of the spectral gap g with 
the system size A^ while its distribution displays Wigner-type level repulsion. When 
the strength of external field goes to zero, and time-reversal invariance is restored, we 
find long range order, logarithmically divergent entanglement entropy, and exponential 
decay of the spectral gap, while the level repulsion disappears. 

We argue that this quantum phase transition is non-conventional from the point 
of view of established models, since in what we shall call non-critical case we still find 
slow power law closing of the spectral gap. 

2. The embedded ensemble of spin Hamiltonians with random two-body 
interactions 

In this section we define the two-spin random ensembles of Hamiltonians for systems 
with N distinguishable spins or qubits with at most two body interactions and describe 
its basic (invariance) properties. If we do not allow all spins to interact, we have to define 
which ones do. The simplest case will be a chain with nearest neighbour interactions, 
but in general we need a graph, whose vertices correspond to spins and whose edges 
correspond to two-body interactions. We proceed to formalize this. 

Let Q = (V, S) be an undirected graph with a finite set of A^ vertices V and a set 
of M edges S C V x V. In addition, let A : V ^ M^ and /i : £^ ^ IR+ be two non- 
negative functions defined on the sets of vertices and edges, respectively. To such a graph 
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we assign a 2^ dimensional Hilbert space of A^ spins or qubits, placed at its vertices 
Ti-g = ®jevC^ = C^ , and a set of A^ Pauli operators a" : Hg -^ Hg, a E {1, 2, 3}, j G V 
satisfying SO{3) commutation relations [o"°, a^] = iSap^a'Jdj^k- We also use the notation 
a, = (a],a|,a3). 

Let A''^-'^^ G M^^^, (j, A;) G £^, be a set of M random real 3x3 matrices, and W^ G M?, 
j G V, be a set of A^ random 3 dimensional real vectors. The TSRE then consists of the 
random Hamiltonians 

The above defined functions on the edges and vertices of the graph are used to 
determine the average strength of the corresponding terms in the Hamiltonian. The 
distribution of random two-body interaction matrices (for short also bond matrices) A^^'^'^ 
and the random external field vectors fe'--'^ shall be uniquely determined by requiring 
the following two conditions: maximum local invariance and maximum independence 
expressed formally as: 

(i) An ensemble of Hamiltonians ([T]) should be invariant with respect to an arbitrary 
local SO {3) transformation, namely 

a'^ = 0,a,, (2) 

where Oj G SO{3)j, j G V, meaning that the choice of local coordinate 
system is arbitrary for each spin/qubit. Obviously, ([2]) preserves the canonical 
commutation relations for the Pauli operators. Then it follows immediately that 
the joint probability distributions of {A^^'''\ W^ should be invariant with respect to 
transformations 

^U,ky ^ oA^^'^'^O', py = 0P\ (3) 

where O, O' are arbitrary independent SO{3) rotations for each j, k. 

(ii) The matrix elements of the tensors A^'J and of the vectors ba should be 
independent random variables. 

Following arguments similar to those presented in [TBj it is straightforward to show that, 
in order to satisfy conditions (i) and (ii) above for pre-determined but general strengths 
of bonds fi{j, k) and external fields A(j), A^'^ and ba should be Gaussian independent 
random variables of zero mean and equal variance, which are uniquely specified in terms 
of the correlators 

(<?e^) =0, 

where (•) denotes an ensemble average. We abbreviate the ensemble defined in this way 
by TSRE(^, fi, A) noting that it depends on the graph and the strength functions fi and 
A. 
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Let us now describe some other elementary properties of TSRE. We have seen 
that each member H of TSRE(^, /x, A) can be parametrized ([1]) by 9M + 3A^ 
independent random parameters. However, we know for the classical ensembles, that a 
parametrization in terms of the structural invariance group, i.e. the invariance group of 
the ensemble and the remaining parameters is very useful. Similarly in the present case 
for an arbitrary set of local SO{3) rotations Oj, the transformation of the parameters 

^0-,fc)' = ojA^^'^^Ok, P^' = Op^\ (4) 

preserves the spectrum and all entanglement properties of H. In fact, the transformation 
dl]) can be considered as a gauge transformation since, composed with the local canonical 
transformation ([2]), it preserves the Hamiltonian H ([T]) exactly. Two Hamiltonians, 
specified by {A^^'''\b^^^ and {A^^'''\b^^^ }, can thus be considered equivalent, and 
the gauge transformation (jl]) defines a natural equivalence relation in TSRE(^,/i, A). 
Therefore, it may be of interest to consider a simplest parametrization of the set of 
equivalence classes, or in other words the orbits of a Hamiltonian on the ensemble under 
the structural invariance transformations. We thus ask, what is a general canonical form 
to which each element H can be brought by gauge transformations (jl]) and how can it be 
parametrized? Let the integer K denote the number of such parameters. The equivalent 
question for the classical ensembles leads to the eigenvalues as canonical parameters, 
since there the structural invariance group is much bigger. 

Let us first consider the simplest connected graph, namely an open one dimensional 
(ID) chain of A^ vertices {1, . . . A^} and M = N — 1 bonds. There it turns out that the 
matrices A^^'^^^' can be simultaneously symmetrized, namely all A^^'^'^^'' = [A^^'^"^^' ]'^, 
by choosing the following gauge transformation 

Oj+i = Rp^, where R^ := ^OJ+i)[t/OJ+i)]^ (5) 

and where 

y^(i,i+i) ^. f/(j,i+i)/;)a,j+i)r^Oj+i)ir^ /g\ 

is a standard canonical singular value decomposition (SVD) of the original bond matrix 
^a,i+i)^ with f/O'J+i), V"(J'J+i) e 50(3) and D^^'^+^'^ diagonal matrices of singular values. 
Since the initial transformation Oi is still free, we can choose it such Oi := U^^''^'^ that the 
first bond matrix is even diagonalized, A^^'^' = D^^'^\ This symmetrization is unique 
provided that singular values of all SVD's (I6l) are non-degenerate which is the case for 
a generic member H. Thus the number of parameters specifying the bond matrices is 
3 + 6(M — 1) = 6M — 3 and in addition to 3A^ external field parameters this gives 
K = 6M — 3 + 3A^ independent parameters. We recover the original set of parameters 
if we add the 3A^ parameters of the group of local rotations. 

Second, we consider the case of a ring graph with A^ vertices and M = N bonds, 
which is obtained from the previous case by specifying the periodic boundary condition 
A^ + 1 = A^. We see that a general H as given in ([T]) can now be symmetrized only if 
the additional topological condition R := R^ ■ ■ ■ R2R1 = 1 is satisfied. Now all but one 
bond matrix can be symmetrized, for example the last one may in addition be multiplied 
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by a topological rotation A^^'^'^' = ^symmetric-^ leading to three additional parameters. 
Therefore in the case of a ring graph we have K = QM + 3A^ independent parameters; 
again adding the 3A^ = 3M parameters of the local rotations we obtain the full set of 
parameters. 

We can now consider the case of a general (connected) graph. From the previous 
examples it is evident that the only crucial additional parameter is the number L of 
primitive cycles, i.e. such cycles which cannot be decomposed into other primitive 
cycles. It is clear that each primitive cycle adds 3 additional topological parameters (or 
one SO{3) topological R matrix) to the 6M — 3 + 3A^ parameters which we would have 
for the case of a tree graph. Therefore we have 

K = 6M + 3N + 3L-3. (7) 

Counting the cycle-contributions and taking the primitivity criterion into account again 
we obtain the the total number of parameters by adding those of the local rotations. 

The above considerations hold for the case of general bond and vertex strength 
functions fi and A. If, however, these functions are degenerate, or even constant, i.e. 
the average interaction strength and field strength do not depend on edges/vertices of 
the graph, then the structural invariance group of the TSRE may be even larger. In 
the latter case this group is obtained as a semi-direct product of a discrete symmetry 
group of the graph Q and the gauge group of local rotations, the latter being the normal 
subgroup. We hope that the considerable invariance properties of the TSRE will prove 
useful in a future analytical treatment of its properties. 

3. Properties of ground states of the TSRE on a ID chain 

In this section we shall only consider the simplest case of a TSRE on a ID chain of 
A^ vertices. In addition, we consider the most symmetric case of constant strength 
functions, say /i(j, j + 1) = 1 and X{j) = A = const. Such an ensemble of random spin 
chain Hamiltonians 

Af-l N 

H=J2ar A(^'^+^)a,^i + A J2 b^'^ ■ ^j (8) 

i=i i=i 

shall be designated as TSRE(A^, A) where we explicitly assume open boundary 
conditions; we shall, however, also consider a ring graph with periodic boundary 
conditions in which case we shall stress this separately. In particular we shall be 
interested in the zero temperature (ground state) properties of TSRE(A^, A). We note 
that due to the large co-dimension of bond strength space the standard perturbative 
renormalization group of decimating the strongest bonds [20] would not work and we 
have at this point to rely on a brute numerical investigation. 

Still, it turns out that most zero temperature properties of TSRE(A^, A) can be 
efficiently simulated using White's density matrix renormalization group (DMRG) finite- 
size algorithm p^ by which spin chains of sizes up to A^ = 80 could at present be 
achieved. One should not forget that all numerical estimates of ensemble average or 
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expectation value of some physical quantity A, which will be designated as {A), require 
averaging over many, say A/'r, realizations from TSRE(A^, A), such that the statistical 



error estimated as a^ ~ y {{A^) — {A) )/Afr is sufficiently small. Due to the lack of 
translational invariance the implementation of the DMRG is non-trivial and was only 
done for a chain with open boundaries. 

We note that for A = 0, any H as defined in ([8]), or even more generally in ([T]), 
commutes with the following anti-unitary time-reversal operation 

f:aj^-a,, H\x=of = f H\x=o, (9) 

so, for odd N, all eigenvalues of H have to be doubly degenerate (Kramer's degeneracy 
[15]). However, as this represents more a technical than conceptual problem for the 
ground state (or ground plane) properties of TSRE(A^, A) we shall in the following 
restrict ourselves to the case of even N. 

3.1. Distribution of the spectral gap 

Let |0),|1), represent the ground state, and the ffist excited state, of H, with 
eigenenergies £"0, and £"1, respectively. It is well known that the crucial quantity which 
determines the rate of relaxation of zero-temperature quantum dynamics is the spectral 
gap g = El- Eq. 

To clarify the difference with respect to TRI, we plot the normalized gap 
distribution dP/dg where g = g/{g) for both cases, together with the theoretical level 
spacing distributions for the Gaussian unitary ensemble (CUE) of random matrices [21] 
and for an uncorrelated Poissonian spectrum (Fig. H]). For the non-TRI case we choose 
A = 1 and observe, to our accuracy, good agreement with the GUE case for two chain 
sizes A^ = 10 and N = 16. Our results suggest that the GUE-like gap distribution, 
exhibiting level repulsion, also holds in the thermodynamic limit N ^ 00; numerical 
results for odd N give the same results. In the case of A = the level repulsion between 
|0) and |1) gradually vanishes as we approach the thermodynamic limit, although no 
conclusive statement can be made about the limiting distribution. For odd A^, the 
ground state is degenerate for A = and the present analysis does not apply. 

3.2. Size scaling of the spectral gap 

Being interested in the thermodynamic limit, it is an important issue to understand how 
g{N) scales with A^. The theory of quantum criticality [IE] states that g remains finite 
in the thermodynamic limit for non-critical systems, and rapidly converges to zero, as 
A^ ^ 00 for critical systems. 

In Figj2]we plot the ensemble averaged spectral gap (g) versus A^ for different values 
of the field strength A. We find a clear indication that in the non-TRI case the spectral 
gap exhibits universal asymptotic power law scaling 

{g)r^N-'^, with r/^ 0.39 ±0.01, (10) 
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Figure 1. Distribution of normalized gap g — g/{g) for a few choices of even N. In 
the non-T-invariant case an agreement with GUE level spacing distribution is obtained 
whereas in the T-invariant case, A = 0, the level repulsion gradually vanishes in the 
thermodynamic limit. 



whereas in the TRI case, A = 0, the asymptotic decay of the gap is faster than a power 
law, perhaps exponential (g) ~ exp(— ^A^), with ^ ^ 0.07 ± 0.02. According to the 
standard theory [16] both cases, A 7^ and A = 0, should be classified as quantum 
critical, however as we shall see later, the case of slow power-law decaying average gap 
( TTOl) has many-features of non-critical systems, such as finite correlation length and finite 
(saturated) entanglement entropy. Therefore we shall, at least for the purposes of the 
present paper, name the case A 7^ as random non- critical (RNC) and the case A = 
as random critical (RC). 

We note that the results for odd and even N are in agreement in the RNC case. 
Also results for the case with periodic boundary conditions up to A^ = 20 show no 
significant difference from the results in Fig. [2] for any A. 



3. 3. Size scaling of the ground state entanglement entropy 

The second characteristics of quantum phase transitions we choose to investigate in the 
TSRE(A^, A), is the entanglement entropy of a symmetric bi-partition of the chain 

S{N,X) = {ti {N/2+i,...N} [(tr|i,...,^/2}|0)(0|)log2(tr{i,...,^/2}|0)(0|)]), (11) 
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Figure 2. Spectral gap scaling with the system size for different values of parameter 
A (see legend). Note an asymptotic scaling oc TV^"'^^, unless A = where faster than 
power law decay of a gap is observed, perhaps asymptotically exponential (see inset 
for a semi- log scale). 



which measures the entanglement in the ground state between two equal halves of 
the chain. It has been suggested in non-random systems [18] that for critical cases 
S oc log2 A^ whereas in non-critical cases S saturates in the thermodynamic limit. 

Indeed, as shown in Figj3l we find for the TSRE(A^, A) that 5* saturates to a constant 
finite S'oo(A) = limAr_^oo S{N, A) for the RNC case A 7^ 0, while in the RC case it grows 
logarithmically 

S{N,X = 0)^clog2N + c', with c^ 0.17 ±0.02 (12) 

We also note an interesting even-odd- A^/2 effect which slowly diminishes as we approach 
the thermodynamic limit. As pointed out in Ref. [22] such an effect is induced by open 
boundary conditions. For periodic boundary conditions the entanglement entropy for 
RNC case is twice as large as in the TSRE(A^, A) with open boundaries. This further 
confirms the conjecture that only short-range correlations around the boundary between 
the two halves contribute to the entanglement. 

We note that our result is essentially different from results for other models, which 
can be obtained by perturbative real space renormalization group [20], for example for 
the disordered critical Heisenberg chain [23], where c = (ln2)/3 and is in general model 
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Figure 3. Entanglement entropy versus chain length N for different values of 
parameter A (see legend). Logarithmic divergence with an estimated asymptotic slope 
S ~ 0.17 log2 N for the critical case A = is indicated with a dashed line. 



dependent 

The fact that the entanglement is reduced in the RNC case with local disorder 
can be explained by chaotic behavior \25\ signalized by the level repulsion in the gap 
distribution. A similar effect can be observed in localization properties where hopping 
of excitations induced by inter-particle interactions is diminished by introducing local 
disorder [26l [27] . This effect of increased localization is useful for successful quantum 
computing and has important consequences for transport properties such as conductivity 

EHl. 



3.4- Correlation functions 

The most direct probe of criticality is perhaps to investigate of long-range order and 
(space) correlation functions. In order to do this we compute the ensemble averaged 
fluctuation of the spin-spin correlation function between two vertices 

C{j,k) = (|(0|aX|0) - (0|a;|0)(0|af|0)p). (13) 

Note that we have to consider average fluctuations of the spin-spin correlation function 
instead of the correlation function itself, since the latter have to vanish due to the local 
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Figure 4. Ensemble fluctuations of the ground state spin-spin correlation function 
C{r) versus distance r, for different values of the parameter A (indicated in the legend). 
Open symbols indicate results for chain length L = 16 and L = 20 whereas closed 
symbols stand for L = 24. 



gauge invariance properties of the TSRE. Because of the local invariance of the TSRE it 
is enough to consider a single type of correlation function, as the RHS of ( TT3l) does not 
depend on indices a,j3 if j j^ k; in fact in numerical computations we average over a,/3 
in order to improve statistics. We consider the Hamiltonian (IHl) with periodic boundary 
conditions. This allows us to average the fluctuation of the correlation over the chain 
and hence C(r) = -^ J2i C{i, i + r). We expect that the results would be qualitatively 
the same for the model with open boundaries and sufficiently large A^, but we obtain 
better statistics in this way. 

Fig. m shows the averaged correlation function fluctuation C{r) for a few choices of 
the control parameter A and of the chain lengths N = 16, 20, 24. The results for chains of 
different lengths coincide for small distances r whereas for larger r flnite-size effects are 
noticeable. For sufficiently large chains it can be conjectured that the fluctuations (or 
effective correlations) asymptotically decay in the RNC case as C{r) x Co2~^/^ with a 
finite correlation length ^ whereas the decay for the RC case is slower than exponential, 
perhaps a power-law, which indicates long-range order ,^ = oo. 
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3.5. Correlation length and the entanglement entropy saturation value 

In Fig. m we observe that the degree of locahzation depends on the control parameter 
A and the correlation function C{r) decays on larger scales as we approach the critical 
point A = which results in a larger correlation length ^. Eventually, the correlation 
length becomes infinite at the critical point A = 0. In Fig. Owe show the dependence 
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Figure 5. The correlation length ^ as a function of the control parameter A. Open 
and full symbols designate correlation length obtained from the best of C{r) with 
cexp(— r/^) on sets r £ {4, . . . , 8} and r £ {5, . . . , 7}, respectively, all for N = 24. The 
inset demonstrates oc — log2 A scaling. Please note that last point at A = 0.1 is rather 
inaccurate due to insufficiently large system size N, so it is not used in a linear fit (full 
line) in the inset, for which other points with A < 1 have been used. 

of the correlation length ^ on the control parameter A as obtained by exponential fit 
of C{r) for a finite size N = 24. Unlike conventional phase transition as in e.g. |17j . 
the correlation length seems to diverge logarithmically as ^(A) ~ — ^olog2A + const 
with ^0 = 0.26 (indicated in the inset of figJS]), even though algebraic scaling cannot be 
entirely excluded with the numerical data that are available at present. 

Large correlation length has strong effect on entanglement. Long range correlations 
demand longer chains for the entanglement entropy to saturate whereas the saturation 
value itself also grows when the critical point is approached. In Fig. [6] we plot the 
entanglement entropy saturation value for various values of parameter A and observe 
similar behavior as for the correlation length. However, unlike the correlation length the 
quantity which diverges logarithmically when A — > is not the entanglement entropy 
S but its exponential. In fact, numerical data for small A show good agreement with 
2ks ^ _ iQg^ ^ _|_ const, where /c = 4.0 (indicated in the inset of figj6]). Note that 2"^ is to 
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Figure 6. Saturated entanglement entropy 6*00 as a function of the control parameter 
A. The inset demonstrates oc log2 log2(A*/A) scaling, where A* = 4, which is indicated 
by a straight line obtained from a best fit to points with A < 1. 



a good approximation proportional to the effective rank, or the Schmidt number of the 
ground state x^ [IHIES], which denotes the number of eigenvalues of the reduced density 
matrix needed to describe the state of the system up to an error e. In fact, the effective 
rank x^, rather than the entanglement entropy, is the decisive indicator of simulability 
by the DMRG method [30] and, we believe, also a relevant quantity in the description 
of a quantum phase transition. 



4. Conclusions 

In the present paper we have defined a two-body random matrix ensemble of independent 
spin Hamiltonians which are invariant under local 5(9(3) transformations and described 
them in a framework of undirected graphs. As the simplest example, we have studied 
a chain with nearest-neighbour interactions in a random external field and observed 
a non-conventional phase transition when the external field is switched off. The 
system is always critical in conventional terminology as it has a vanishing gap in the 
thermodynamic limit in all cases studied. Yet we have shown that, in the presence of 
a random external field breaking time-reversal invariance, the locally disordered system 
has many properties of non-critical systems such as finite correlation length and finite 
bipartite entanglement entropy in the thermodynamic limit, whereas the gap decay 
obeys a universal power law dependence. The transition towards the critical point with 
vanishing of the external field exhibits logarithmic divergence for the correlation length 
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and the effective ranlc of the ground state. We have no explanation for the logarithmic 
behavior in the quantum phase transition. 

The model proposed is much richer than the example discussed. Thus we expect, 
that higher connectivity of the graph will yield very different results, but even an 
exploration of high temperature behaviour for the chain seems very worthwhile. In view 
of the large structural invariance group of the ensemble in the case of site independent 
average coupling and external fields we hope, that some analytic results can be obtained 
for this ensemble. 
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